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We have made current-voltage (IV) measurements of artificially layered high-T c thin-film bridges. 
Scanning SQUID microscopy of these films provides values for the Pearl lengths A that exceed 
the bridge width, and shows that the current distributions are uniform across the bridges. At high 
temperatures and high currents the voltages follow the power law V oc I n , with n — $o/87r 2 AfcsT-|-l, 
and at high temperatures and low-currents the resistance is exponential in temperature, in good 
agreement with the predictions for thermally activated vortex motion. At low temperatures, the 
IV's are better fit by InV linear in I~ 2 . This is expected if the low temperature dissipation is 
dominated by quantum tunneling of Pearl vortices. 

PACS numbers: 74.50.+r, 74.78.Fk, 74.78.-w 



Soon after the seminal work of Caldeira and Leggett 
0, the possibility of quantum tunneling of vortices in 
thin films was considered by Glazman and Fogel 2] . This 
idea was studied extensively in theory and experiment in 
relation to finite flux-creep rates in hard type-II super- 
conductors for T — > 0. For thermally activated creep, 
the rates should go to zero in this limit, see review |3| 
and references therein. Within creep models, the barri- 
ers through which the vortices tunnel are due to material 
disorder and are treated statistically. 

Narrow current carrying thin-film bridges offer a 
unique opportunity to study both thermal activation and 
quantum tunneling through the well defined potential 
barriers that are tunable by current and temperature. 
Within the London approach, the barrier shape is well 
known for bridges narrow relative to the Pearl length 
A = 2 A 2 /d, where A is the penetration depth and d is the 
film thickness. 

To realize these conditions we use ultra-thin high-T c 
films with extremely large Pearl lengths A ~ 100 /im. 
The critical temperature of these films T c ~ 40 K pro- 
vides a relatively easy T-domain to work in, and the data 
reported are quite robust. 

When current / flows through the bridge at high tem- 
peratures, vortices are thermally activated at the strip 
edges and pushed in, causing dissipation and a non-zero 
voltage V. We show that these processes should lead to 
a power-law dependence V oc I n with the exponent n 
determined by the film parameter A and temperature: 

(</>o is the flux quantum and k is the Boltzmann constant). 
This is confirmed by the high-temperature IV's and by 
the independently measured Pearl length. 



However for low T's, the IV's show a slower current de- 
pendence than required by thermal activation. Instead,a 
quantum tunneling model utilizing the ideas of Refs.[l|,|2j, 
along with our knowledge of the barrier shape, provide a 
good representation of the data. 

In the present work, ultrathin sandwiches of 
[Bao. 9 Ndo.iCu0 2+;c ] 5 /[CaCu0 2 ] 2 /[Bao.9Nd . 1 Cu0 2+:c ]5 
(CBCO 5/2/5) were used. They consist of one super- 
conducting infinite layer block (two CaCu0 2 unit cells), 
sandwiched between two charge reservoirs (5 Ba-based 
unit cells). The structure was grown on (001) SrTi03 
substrates, with nominally zero miscut angle, by using 
a focussed KrF excimer pulsed laser source |4j, |5J . The 
films were photolithographically patterned into bridges. 
In this paper we present measurements from a "wide" 
bridge approximately 85/im across. After measurement, 
this bridge was narrowed by a factor of 2 and then 
remeasured. The patterning of such a thin layer (5 nm 
with a superconducting cell only 1 nm thick) protected 
by a 100 nm amorphous cap layer required specially 
careful photolithography and milling. The use of a 
cap layer deposited in situ required pre-deposited Au 
contacts (details will be given elsewhere). Four terminal 
measurements of the current- voltage characteristics were 
realized down to 300 mK in an Oxford Heliox probe 
suitably modified by superconducting wiring to minimize 
heating at elevated currents. The bridges were imaged 
with a scanning SQUID microscope (SSM) H to 
determine current uniformity and the film Pearl length. 

Figure nji shows an SSM image of the wide bridge af- 
ter cooling in a 22 mG field. The bridge outline is visible 
along with trapped Pearl vortices HjEJ. Fig. \I]p shows 
a cross-section through one Pearl vortex, after recooling 
in a smaller field to reduce the density of vortices. Fit- 
ting such cross-sections results in low temperature Pearl 
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FIG. 1: (Color online) (a) SSM image of the wide 5/2/5 
CBCO bridge after cooling in a 22 mG field at 5 K with an 
8 /im square pickup loop, (b) Cross-section through an iso- 
lated Pearl vortex after cooling in a 1 mG field. The line 
fit is obtained using the field distribution h z (x,z) given in 
Ref.H with an effective pickup loop height 2 = 5 fj,m and 
A = 400 ^m. (c) SSM image with a 104 Hz, 100 /iA r.m.s. a.c. 
current passing through the bridge at T =20 K. (d) Cross- 
section through the data along the dashed line in (c); the solid 
line is calculated for uniform current through the bridge, with 
2 = 5 /jm. 



lengths of 400±50 /im for this bridge. Fits to scanning 
susceptibility measurements of the same bridge gave A = 
200±20/im (see Ref.0 for details of these procedures). 
This difference may be due to the finite lateral dimensions 
of the film, which are not accounted for in the modelling. 
The length A is therefore longer than the bridge width, 
and we expect the current through the bridge to be uni- 
form, as confirmed by the SSM image and fit of Fig.^,d. 
The field h z (normal to the film) is maximum near the 
edges x = 0, W where vortices or antivortices nucleate, 
penetrate the strip by being pushed in by the Lorentz 
force 4>qI/cW, and annihilate near the strip middle. We 
focus on the barrier for vortices; the one for antivortices 
is obtained by replacing x with W — x. 

For narrow strips of width W <C A with no current, 
the energy of a vortex at a position < x < W is : 
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£ is the coherence length. The vortex energy in the pres- 
ence of a uniform current is 



U(x) = e(x) 
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(3) 



The current causes suppression of the barrier maximum 
U m and pushes the maximum position x rn from the mid- 
dle (at 1 = 0) toward the edge x = 0. We obtain after 
simple algebra: x m = (W/ir) tan" 1 (/ /^) and 
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For A « 400 jum, I k2x 10~ 3 mA at low temperatures 
and decreases on warming. For the majority of our data 
7/ Jo > 1 and Eq. |g} simplifies to: 
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where e = 2.718. The logarithmic dependence of the 
barrier height on the current results in power-law IV's, 
the dependence encountered in models of vortices moving 
through random potentials in disordered materials Q . In 
our case, it is caused by a well-defined barrier. 

If vortices cross the strip via a thermally activated pro- 
cess, a voltage V oc exp(—U m /kT) is generated, i.e., 



V oc exp [ — -=ln -y 
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where Id is of the order of the bridge depairing current. 
Thus, the model provides not only the power-law per se, 
but gives the exponent in terms of the film parameters. 
We note also that our model differs from that used in 
RefJUl to interpret similar data. 

As is usually the case with thermal activation, it is 
difficult to fix the "attempt frequency" in the probability 
exp(— U m JkT). Still, the number of vortices penetrating 
the strip must be proportional to the field value at the 
edges, i.e., to the current I. Adding a prefactor const-/ 
in Eq. Q , we obtain the exponent n = m + 1 of Eq. |QJ 
which gives the correct limit for T — > T c . 

The upper panel of Fig. |2 shows IV measurements for 
the wide bridge on a log-log scale. All curves are nearly 
straight (except noisy low- voltage parts), i.e., V oc J™ 
at high currents. The exponent n extracted from power- 
law fits to the data is shown in the lower panel, for both 
the wide and narrow bridges. The solid line is n(T) of 
Eq. Q), in which the two-fluid A = 2\^/d(l~t A ) has been 
used for simplicity; t = T/T c . Fitting the high-T part of 
the data for the wide bridge we obtain T c w 38 K and 
A(0) « 320 /im in agreement with A(0) w 200 - 400 /um 
independently measured by SSM. Thus, the thermal ac- 
tivation model works well above ~ 15 K, where it gives 
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FIG. 2: (Color online) The upper panel: V(I) on a log-log 
scale at temperatures indicated in the legend. Lower panel: 
The exponent n extracted from V oc I n for the wide bridge 
(circles) and the narrow bridge (triangles). The solid line is 
n(T) of Eq. £Q| with with the best-fit parameters T c — 37.6 K 
and A(0) = 317 ^m. 



correct values of the exponent n(T). Although T c of 
the narrow bridge was slightly lower after repatterning, 
the high temperature exponents were quite close, which 
means that the Pearl length was nearly the same in both 
bridges. 
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FIG. 3: (Color online) R(T) for 30 < T < 36 K and fit to an 
exponential dependence. 

Another quantity which can be extracted from the IV's 
is the low-voltage differential resistance, defined as R — 
(dV/dI)i-*o. At low temperatures, the voltages at low 
currents are exceedingly small and cannot be extracted 
from the noise. However, reliable determination of R 



is possible above roughly 30 K. Fig-El shows that in this 
domain R oc exp(T/Xo). 

To understand this result, we turn again to the volt- 
age expression Q for thermally activated vortex mo- 
tion. As mentioned above, the prefactor in this expres- 
sion is proportional to the current because the number 
of vortices attempting to overcome the barrier is propor- 
tional to the edge fields, i.e., to I. Therefore, calculating 
R = (dV/dI)i^o, we can leave the current dependence 
only in the prefactor and set I = in exp[— U m (I)/kT]. 
We then have: 



2W _ 4<£§(1 -t) 2W 
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where A(0) is the low temperature Pearl length, and the 
two-fluid approximation A = A(0)/(1 - 1 4 ) w A(0)/4(l - 
t) is again used. Then, the voltage and the resistance 
near T c should be proportional to 
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where To depends weakly on T via ln£. Note that our 
model does not allow us to literally take the limit T — > T c 
because the energy expression (|2J is valid only for W ^> £. 
The data of Fig.[3]correspond to To w 0.78 K. Expression 
(|12fl yields a smaller value To « 0.3 K. 

We stress that R oc exp(T/To) is obtained here for the 
thermally activated process, so that this dependence per 
se is not necessarily a manifestation of quantum tunnel- 
ing (compare with (l2^h 

At low temperatures, the thermal activation model 
fails. The fit values for n(T) in Fig. [21 do not diverge 
as thermal activation requires [l3j |. Hence, we turn to 
the possibility of quantum tunneling. According to Ref. 
2, the tunneling probability for overdamped processes is 
proportional to exp (—jijxl /m where 7 is a geometric 
factor related to the barrier shape, 77 is the drag coeffi- 
cient, and Xb is the barrier width. Unlike thermal activa- 
tion, for which only the barrier height is relevant, for the 
tunneling the barrier width dominates the probability. 

Since the potential U = at the position of vortex 
entry, the width Xb is a root of U (x) = 0, or of 
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As argued above, I/Iq 3> 1; therefore irx/W must be 
small. Since the log is a slow function, we have 
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where Co is a constant of the order one. Hence, the 
barrier width shrinks with increasing current. One can 
show by solving numerically Eq. (| 1 31) that for I/I ~ 10 3 
(the range of our data) C ~ 5. 

We then expect the voltage for high currents to behave 
according to 

^V = Cl -C 2 ^f§. (15) 

Here, C\ is a constant related to the attempt frequency, 
whereas C 2 = Cqj/tt 2 ~ 2.57. 

In Fig.0J we plot \ogV for a few T's versus I~ 2 for both 
wide and narrow bridges. Clearly, the straight parts of 
the IV's for low temperatures and high currents are in 
agreement with Eg. (|15l) . With increasing T, the IV's 
curve away from a straight line much faster, signalling 
deviations from predominantly quantum tunneling. We 
fit the data for 0.8 K of the wide bridge and for 1 K of the 
narrow one to V = Vb exp(— a/I 2 ) to obtain Vb=0.43V, 
a = 1.0 x 10~ 4 A 2 for the first and V o =0.97 V, a = 0.14 x 
10 _4 A 2 for the second. 
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FIG. 4: (Color online) The voltage V versus l/I 2 for tem- 
peratures indicated in the legend. The left group of IV's is 
for the wide bridge with W = 85 urn; the right group is for 
IF = 42.5 /mi. 

Going back to Eq. (|15l) for the quantum tunneling, we 
note that the coefficient of I~ 2 is proportional to W 2 . 
Therefore, as an extra check of applicability of this equa- 
tion, one can compare a's of the two bridges. Given the 
approximate character of our model, the reduction of a 
by a factor of 7 in the bridge W/2 as compared with the 
expected reduction by a factor of 4, shows that quantum 
tunneling is consistent with our data. 

Moreover, the pre-exponential factor Vo proportional 
to the attempt frequency should go as 1/W because - 
as magnetostatics shows - the edge field is oc 1/W. The 
ratio of Vo's for the two bridges is 0.97/0.43=2.3 instead 
of the expected factor of 2. This is yet another indication 
that the low-T data are consistent with the tunneling 
model. 

Further quantitative tests of the quantum tunneling 
model would be possible if the drag coefficient n and the 



barrier shape factor 7 were well known. Unfortunately, 
this is not the case. The factor 7 = ir/2 was found for a 
barrier U(x) as a cubic parabola |l5|. quite different from 
our case. For I — 2.5/iA and W=85/im, we estimate 
777 ~ 10~ 16 CGS. If one takes 7 ~ 1, this gives the drag 
coefficient about two orders of magnitude smaller than 
the Bardeen-Stephen value ij = (fy^d / '2tt^ 2 c 2 p n . Kramer- 
Pesh contraction 0] of the vortex core at low T's might 
be a cause of the n suppression. 

During this work, we have also searched for a 
Berezunskii-Kosterlitz-Thouless transition in the domain 
of the power-law IV's. We, however, did not detect any 
discontinuities in the T dependence of n. We believe this 
is because in bridges with W < A, interaction of vortices 
is short-range with a cutoff at distances W. 

To summarize, at high temperatures we observe the 
power-law IV's well described by thermally activated mo- 
tion of vortices through the thin-film strips narrow on 
the scale of Pearl A for which the shape of the potential 
barrier is well established. The predicted dependence 
of the exponent n on A(T) is in good agreement with 
the data. The low temperature data are consistent with 
the overdamped quantum tunneling of Pearl vortices, the 
size of which in our films is indeed macroscopic since 
A ^ 0.1 mm. 
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